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Abstract. The process of un-reduction, a sort of reversal of reduction 
by the Lie group symmetries of a variational problem, is explored in 
the setting of field theories. This process is applied to the problem of 
curve matching in the plane, when the curves depend on more than 
one independent variable. This situation occurs in a variety of instances 
such as matching of surfaces or comparison of evolution between species. 
A discussion of the appropriate Lagrangian involved in the variational 
principle is given, as well as some initial numerical investigations. 
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1 Introduction 

The idea of un-reduction was introduced in for the purpose of using reparametri- 
sation by the action of the group G = Diff(S'^) to improve resolution of selected 
features in dynamics and optimal control problems, particularly in matching 
closed planar curves, whose configuration space Q is the space of embeddings 
Q = Emb(5'^,M^), by redistributing grid points in along the curve. The 
un-reduction process used in [9] was based on reconstruction, the inverse of La¬ 
grangian reduction by symmetry m, which relates the solutions on TQjG to 
solutions on TQ. This version of the un-reduction process was applied in the 
outer metric setting in m- In this setting, the deformation of the shape of 
a curve in Emb(5'^,M^) was applied to the embedding space, K^, completely 
independently of any reparametrisation of the embedded space, S^. 

In contrast, the un-reduction technique introduced in [1] seeks a family of 
equations, called the un-reduced equations, on TQ, whose solution projects onto 
those of a set of Euler-Lagrange equations on T{Q/G). Thus, the un-reduction 
process used in [4] is distinct from the reconstruction process used in |9I8) . More¬ 
over, the un-reduction approach introduced in [1] raised an important issue for 
numerical applications in curve matching by optimal control, since it intertwined 
the reduction and reconstruction processes. Namely, the geodesic distance be¬ 
tween two curves should be independent of their parametrisations. In particular. 


2 


Arnaudon et al. 


the measure of the deformation of the shape of a curve should be independent 
of its parametrisation. 

Here we address this issue for curve matching from the viewpoint of a re¬ 
formulation of the un-reduction scheme of [3], which gives a framework in the 
setting of inner metric discussed in [3] , in which the shape deformation is applied 
to the embedded space, rather than the embedding space. 

In the process of addressing this issue, we will also generalise the un-reduction 
scheme by formulating it as a covariant space-time field theory. This generalisa¬ 
tion gives us the freedom to introduce additional independent ‘time’ or ‘space’ 
variables for the purpose of coordinating comparisons among shapes. Introduc¬ 
ing additional independent variables allows more flexibility in making shape 
comparisons than, for instance, the time warp approach of |10] . which does not 
increase the number of independent variables. For example, one could imagine 
making comparisons of cylindrical surfaces by assembling closed curves resolved 
on two-dimensional slices. In this case, the additional space variable would be 
transverse to the slices, and one would make comparisons of surfaces as single 
entities, rather than comparing the evolution in time alone of each slice indepen¬ 
dently. In addition, the covariant field theory generalisation of the un-reduction 
framework in the inner metric setting could lead to a variety of other applica¬ 
tions, a few of which are mentioned in Section]^ 

This work may be summarized as an extension of [3] in the following three 
directions: 


( 1 ) 


( 2 ) 


We promote the un-reduction formulation of [3] in classical mechanics to a 
covariant field theory by following the same reasoning. Namely, we derive the 
Lagrange-Poincare reduction of the shape space (Section 
it by including an independent parametrisation (Section 


3.1) and un-reduce 


3.2). 


Instead of the curvature weighted metric used in [3] , we implement Sobolev 
metrics, which avoid the issue of arbitrary small geodesic distances (Section 


^ 1 . 

(3) We finish by illustrating this approach and assessing its validity with a few 
numerical experiments in the classical mechanical setting (Section |4.1[). 


A complete exposition of covariant un-reduction containing the proofs and ap¬ 
plications in other areas, such as in theoretical physics, can be found in [T]. 

The main topic of the present paper is the first point listed above. We shall 
focus our discussion on the covariant Lagrange-Poincare (LP) reduction by sym¬ 
metry in the context of curve matching. The symmetry will be G = Diff(S'^), the 
diffeomorphism group which acts on the configuration space of planar curves, 
Q — Emb(S'^,M^), as a reparametrisation of S^. LP reduction is a general 
method to deal with noncanonical reductions, in which the configuration space 
is not the symmetry group. LP reduction allows the explicit derivation of the dy¬ 
namical equations on the quotient space, Q/G, which, in our case, is the space of 
shapes, Emb(S'^,K^)/Diff(S'^). The field theoretic version of LP reduction that 
we will use here was developed in mu, based on the classical reduction theory 
introduced in [7] . We will use a simplified version of this theory applied directly 
to curve matching. 






Covariant un-reduction for curve matching 


3 


The second improvement relative to [4] on the list above involves the Rie- 
mannian metric used to derive the un-reduction equations. As recently pointed 
out in |3I2) , use of Sobolev inner metrics avoids the problem of having arbitrar¬ 
ily small geodesic distances between two curves. In addition, the un-reduction 
equations take a simpler form with the Sobolev metric than with the curvature 
weighted metrics used in [i]. 

Finally, we tested the un-reduction approach in a few numerical experi¬ 
ments where we considered an initial value problem with reparametrisation. 
Even though we simply chose the forward Euler method for this initial value 
problem, without any further modifications, it still converged to the expected 
solution. The success of these simple numerical simulations motivates us to go 
further in future work to consider boundary value problems in the full field 
theoretical framework as explained in Section [4.2| 

2 The geometry of curve matching 

We start by recalling some basic facts about the geometry of curve matching 
that we will use throughout the text. We refer to earlier works such as |14 |3|2|4j 
for more details. 

2.1 Reduction structure. 

Let Q = Emb~''(S'^,K^) be the manifold of positive oriented embeddings from 
to K^. Elements in Q are maps c(0) € for 0 € and elements in the 
tangent space T^Q are pairs (c,v) with c S Emb'''(S'^,and u € 
a parametrized vector field along the curve c, thus TQ = Q x 

We then consider the group G = Diff'''(S'^) of orientation preserving diffeo- 
morphisms of and its Lie algebra g = X(S'^). The group G acts on the right 
in Emb''’(S'^, as reparametrisation of curves c and the reduced space is the 

space of shapes in 


Q _ Emb+(S'\K2) 
G “ Diff+(S'i) 


The reduction of the phase space TQ hy the group of symmetry G is a compli¬ 
cated space but can be decomposed via the introduction of a principal connection 
A : TQ —>■ Q as follows 


G 


A 


■> TA © 0 = T 



Q X 0 
G 


( 2 ) 


The space FA is then the tangent space to the space of shapes, and the adjoint 
bundle 0 will encode the parametrisation velocity of the curve. This space seems 
rather abstract but it corresponds to having a Lie algebra attached to each point 
of the base space A, which means that for each shape we have the freedom to 
attach an arbitrary parametrisation velocity. The construction of the connection 
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is straightforward for curve matching. Namely, given the velocity u S T^Q, we 
consider its tangent and normal decomposition, 

u{9) = v{0)t{O) + h{9)n{0 ), (3) 

where (t, n) is the orthonormal Frenet frame along c and v(9), h{9) G X(S^) are 
scalar functions along the curve, parametrised by 0 S We clearly have that 
v(9)t(9) is a vector tangent to the orbits of G = Diff'''(S'^). This decomposi¬ 
tion defines the principal connection A : TQ —>■ g, which, when applied to u, 
gives A{u) = V, the reparametrisation velocity. The horizontal part of u is then 
h{9)n{9) and we have a decomposition TQ = HQ 0 VQ. We will also need the 
curvature of A, defined as B := d-^A = dA + [-4,^], but its exact form will not 
be needed here. In fact, we shall skip any technical details which are not directly 
useful for the present work, and refer the interested reader to mm for the full 
discussion of this construction. 


2.2 Field theoretical structure. 

We can now extend this reduction structure by promoting the classical system to 
a covariant field theory, see mu . In order to do this, we consider an open domain 
TV C K X K endowed with the Euclidean metric, the associated coordinates (t, x) 
and a volume form v = dtdx. For simplicity we only consider a two dimensional 
space-time manifold N, but more dimensions can be added in a straightforward 
way. The tangent space TQ is then promoted to the jet bundle in order to capture 
the space-time direction used to compute a tangent vector. In this simple setting, 
the jet bundle has a simple geometric meaning, given by J^{N, Q) ~ T*N ®TQ 
and a generic element will be written 

j^c = ct{9){t,x)dt + Cx{9){t,x)dx, (4) 

that is, Ct,Cx G TQ are the derivatives of a map c : N ^ Q along t and x 
respectively. The generalisation of the time derivative is the divergence operator, 
given in this simple case by div = ^Ct + -^Cx- 

2.3 Riemannian metrics 

The choice of a convenient Riemannian metric on Q which is invariant with 
respect to the action of G is an interesting topic. See, for example, [5] and [3] 
and the references therein for more discussion. In our case, invariance under 
reparametrisation is achieved by considering arclength integrations and deriva¬ 
tives ds = \ce\d9 and Dg = \^\dg. However, the main difficulty lies in the 
geometrical properties of the metric, as we deal with an infinite dimensional 
space. The natural metric 



9l^{u,v) 


{u{9),v{9))ds, 


(5) 
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with u,v G TcQ such that (•, •) is a dot product in is not very useful as it can 
lead to arbitrarily small geodesic distances in both Q and Q/G, see [H], The 
problem can be overcome in the shape space Q/G by the metric 

9 k(u,v)= [ {1 +AK{6f‘){u{e),v{B))ds, ( 6 ) 

JSi 

where ^ > 0 and k, denotes the Frenet curvature of the curve c, defined as 

K := (DeDgc) ■ J{Dgc) = Dgt • n. (7) 

In fact, the weighted metric in ([^ can still have arbitrarily small geodesic dis¬ 
tance in Q along the fibres of the fibration Q —>■ Q/G. A metric with a well 
defined Riemannian distance in both Q and Q/G may be obtained by adding 
higher order derivatives of u and u in a Sobolev-type expression as 

(u, u) = / {{u{9),v{9)) + A‘^{Dgu{6),D0v{9))) ds. (8) 

We can collect these three cases (as well as others, see [5]) by using 

gv{u,v)=f {u{9),Vv{9)), (9) 

Js^ 

for a convenient choice of a G-invariant self-adjoint pseudo-differential operator 
V which can depend on the curve and its derivatives. In particular, the operator 
for Q is 7^ = 1 -I- Ak^ and for Q we have 7^ = 1 — A^Dg. One additional 
advantage of the operator associated to Q is that it does not depend on the 
curve, whereas the operator for depends on the curvature of the curve where 
it is evaluated. This represents a great simplification in the expression of the 
un-reduced equations. 

3 Reduction and un-reduction 

We are now ready to perform reduction by symmetry from the space of embed¬ 
dings to the shape space using the covariant Lagrange-Poincare reduction. 

3.1 Lagrange-Poincare reduction 

Let’s begin by recalling the original problem of curve matching. The matching 
problem is a boundary value problem in A = Q/G with Lagrangian I : TS ^ M.. 
Hamilton’s principle states that the Euler-Lagrange equations associated with 
this Lagrangian yield the solution which minimises the action functional given 
by 1. In practice, the matching problem is solved by using a shooting method for 
determining the initial momentum such that the curve at the final time matches 
the target curve within some specified tolerance. Instead, we will start in the 
larger space Q, where the numerical experiments can be done easily and, as a 
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first step, reduce this system such that we recover the Euler-Lagrange equations 
on S. 

We project the variational principle defined for L from J^{N,Q) to its quo¬ 
tient J^[N, Q)/G = J^{N, S) 0 {T*Ni^g), where § = (Q x 0 )/G, as in equation 
(|^. Critical solutions are maps a : N ^ T*N which, moreover, project to 
maps p:N^E = Q/G as p = TTgOa according to the diagram 

T*NiEs 

y ( 10 ) 

N E 

where TTg : T*N (Si q ^ E is the projection of the adjoint bundle neglecting 
the T*N factor. The free variations of the initial problem provide a family of 
constrained variations that define a new type of variational equations, called 
Lagrange-Poincare equations, [5], m- The next theorem gives the Lagrange- 
Poincare reduction which includes forces F : T*N ^TQ —>■ TQ. Before stating 
this theorem without proof, we will make another important assumption which is 
satisfied by most of the Lagrangians used in curve matching. Namely, we assume 
our Lagrangian decomposes as a sum of two Lagrangians taking values from the 
vertical and horizontal space that will be denoted L — + L" and i + F. 

Theorem 1 (Covariant Lagrange-Poincare reduction with forces). Given 
a map c : N ^ Q, let a : N ^ T*N 0 g 6 e defined as 

CT(a:)(a;) = [s(a;),^(Ts • (w))]g , (11) 

with uj S TxN,x S N and where [-Jg stands for the quotient by G; p : N ^ E, 
p{x) = [s(x)]g = TTg o (t. With the previous definitions, the following points are 
equivalent 

(1) s is a critical mapping of the variational prineiple 

^ [ L{s,j^s)v+ [ {F{sJ^s),Ss)v = 0 (12) 

Jn J n 

with free variations Ss. 

(2) The Euler-Lagrange form of L satisfies the relation 

£C{Lv)ifs) = F, 

where £C is the Euler-Lagrange operator acting on the seeond jet bundle 
(second order field theoretical tangent space), whieh gives the usual Euler- 
Lagrange equations. 

(3) a : N ^ T*N^ q is a critical mapping of the variational principle 

S [ £{j^p,a)v+ [ {f^(j^p,a),5p)w F [ p,a),ri)v = 0, 

Jn Jn Jn 

for variations of the form Sa = Er] — [a,r]] + B{Sp,Tp) G 0 , where 6p G TpE 
is a free variation of p and rj is a free section 0/0 —> 17. 
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(4) cr satisfies the Lagrange-Poincare eguations, written if L and I decomposes 
according to the vertical/horizontal decomposition 


div 


St 

Sfp 


5p ■’ 5p 

<yr ,*(5r , 


SP’ 

Sa 


,B{Tp,-) 


(13) 


One recognises left hand side of the first equation in (131 as an Euler-Lagrange 


equation and the second one as an Euler-Poincare equation. The right hand side 
of both equations are either forces, or coupling with between them. The solu¬ 


tion a of the Euler-Poincare equation in (13) will influence the Euler-Lagrange 


equation via the term involving the curvature of the connection A. An addi¬ 
tional coupling arises because a is in the adjoint bundle and therefore depends 
implicitly on the base curve in Q. 


3.2 Un-reduction 


The particular form of the equations in (13), based on the decomposition of the 


Lagrangian and the inclusion of the force term will allow us to decouple these 
equations in the sense that the right hand side of the EL equation will vanish; 
so the feedback of the EP equation to the EL equation will disappear. Before 
stating the un-reduction theorem we must recall that the canonical momentum 
map J : T*Q —i' g* for the natural lift action of G on T*Q, is defined by 

(J(a,),^)gxg* = {aqAQ)TQxTQ- , 

where Oq G T*Q, € g, and fg € TQ is the infinitesimal transformation of 
the action of G on Q at the point q € Q. The map J extends to a map J : 
TN ®T*Q TN ® g*, trivially in the factor TN. 

We can finally state the covariant un-reduction theorem. We refer to [T] for 
the proof and more details about this theorem. 

Theorem 2. We consider a G-equivariant force F : J^{N, Q) —>■ T*Q such that 
F"" = p" o F is arbitrary and F^ = p^ o F is given by the condition 


Vf 


r = -^ + {^MTp,-) 


Sp 


St 


5a 


(14) 


for its projection f^ : J^{N,S) x (T*A^ 0 g) —> T*S. Then, the variational 
equations of the problem defined by L and F read 


SC{L^){j^s) = 0 


(15) 


where A* : g* V*Q is the dual of the connection form. Finally, critical so¬ 
lutions s : fV —)■ Q o/ ( 15 ) project to critical solutions p = [s]g of the Euler- 


Lagrange equations £C{l){j^p) = 0. 
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Remark 1. For TV = M, v = di, that is, in the case of classical mechanics, we 
have div = d/dt and we recover the results and equations of [4]. 


The first equation in (151 is the usual Euler-Lagrange equation for the horizontal 
Lagrangian needed for solving the matching problem. Regarding the interpre¬ 
tation of the second equation, the definition of J above shows that 3{5L'"/5j^s) 
is a covariant momentum map, so that div J((5L’'/(5j^s) is the expression of a 
conservation law with respect to the group of symmetries. If one set F'" = 0, the 
conservation law is complete. However, sometimes it may be interesting to keep 
this vertical force, as it might be used to externally control the dynamics along 
the vertical space; that is, the reparametrisation. 


3.3 Un-reduction with Sobolev metric 

We consider the Diff^(iS'^)-invariant Lagrangian L : J^{N, Q) ~ T*N®TQ M 
which can be decomposed as L = + L'" with respect to the connection A as 


^ / {ht'Pht + K'Phx)ds, L'’{j^c) = ^[ {vtVvt + v^Vv^) ds 


where 


Ct = Vft + htii and Cx = Vxt + hxn. 


The un-reduction equations (151 are then computed in Proposition below in 
the case when V is independent of the curve, that is, Sobolev metrics. 


Proposition 1. The un-reduced equations (15) for the hi-dimensional problem 
of planar simple curves defined by the Lagrangian defined above and the metric 
§ with V being the Sobolev operator, read 


dxVhx + dtVht = DeihxVvx htVvt) - kH 
dxVvx + dtVvt = F" 


(16) 


for any choice of vertical force F'", where 


id — - {hxVhx + htVht) ■ 


(17) 


In (16), the function FI is the shape kinetic energy density and n is the Frenet 
curvature of the curve c, defined in 0 - This term can be interpreted as a penalty 
term in deforming curved regions. The sign of this term would depend on the 
concavity or convexity of the curve at this point, and thus this force would try to 


prevent the curve to be deformed too fast in these regions. Equation (16) shows 


that the dynamics in (x, t) is governed by the coupling between ht and Vt required 
for the shape deformation to be independent of the reparametrisation. In fact, it 
also contains a derivative with respect to 0 and vertical vectors. This is the only 
term which couples with the vertical equation, and it also gives the corrections 
necessary for the curve deformation to be independent of the reparametrisation. 
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Remark 2. In the classical mechanics setting, the un-reduction equations with 
Sobolev metric would be very similar to (16), but without the x-dependent terms. 

Owing to the 


They will be the equations used for the numerics in section |4.1| 
simplicity of the Sobolev metric compared to the curvature weighted metric, the 
derivation of this equation is directly done from the un-reduction equations, not 
from the variational principle, as in [i]. We refer to [T] for the details of this 
calculation. 


4 Applications 

Before discussing the possible applications of the covariant un-reduction scheme 
in curve matching, we shall present a short numerical study of the un-reduction 
in classical mechanics using the Sobolev metric . 

4.1 Numerical validation 

In order to test and illustrate the un-reduction scheme, we performed a few 
simple numerical experiments. We restricted ourselves to the classical case, al¬ 
ready done theoretically in but with the norm instead of the curvature 
weighted norm. The only effect of the Sobolev norm that which interested us is 
that it regularises the curve deformation and prevents large bending of the curve, 
smaller that the scale given by A. A deeper analysis of the effect of the Sobolev 
norm in the matching process is not the aim of this paper but is important for 
applications. Our main goal here was to check the decoupling between the shape 
and the reparametrisation dynamics for a simple initial value problem. 

Our numerical scheme used the Euler explicit scheme in time and 2nd order 
centred finite difference approximation for Dg in order to have a symmetric 
space discretisation in 6* o —9. The application of the Sobolev operator V = 
1 - was done in Fourier space with A = 0.3. Our initial condition was 

a circle and the initial horizontal velocity was a bump function, so the curve 
deformed as in Fig. [l] The curve was discretised with 100 points and we used 
a set of time steps {dt = 0.04,0.02,0.01,0.001) to study the convergence of the 
scheme, and especially the decoupling between the reparametrisation and the 
shape deformation. In order to do this, we ran two initial value problems where 
one of them also had a vertical initial velocity, taken to be constant such that 
the parametrisation would rotate during the evolution of the curve. We then 
computed the distance between the two curves using the methods of currents 
|12j at each time to make a parametrisation-free comparison of the shapes of the 
curve. The results are displayed in Fig.j^together with the distance between the 
curves as a function of time. Even with the simple numerics we used (the Euler 
scheme and finite difference), the un-reduction feature was verified. The example 
we studied is simple and did not require high resolutions and robustness tests 
as for more realistic shapes. Further numerical studies using the un-reduction 
scheme would thus include improvements of the current implementation and 
a shooting method in order to solve the correct matching, or boundary value 
problem. 
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X 

(a) dt = 0.04 


X 

(b) dt = 0.02 




X 

(c) dt = 0.01 (d) Distance between fixed and 

reparametrised curves with the method of 
currents 

Fig. 1. This figure illustrates the effect of the time step of 5''^-reparameterisation on the 
quality of the deformation of a curve using the explicit Euler scheme in the un-reduction 
approach. The blue curve has horizontal initial velocity ht for shape deformation, while 
the green curve has an additional constant vertical initial velocity vt, for reparameter- 
isation. The blue and green dots begin at the same initial point, but then the green 
one shifts along the green curve as the reparameterisation proceeds. The black dashed 
curve is the initial condition of the simulation. Upon decreasing the time steps, the 
coupling between the vertical and horizontal dynamics decreases and the quality of the 
deformation improves, even when using an explicit Euler scheme. On the bottom right 
panel, we plotted the distance between the two curves, as a function different time 
resolutions. The distance is computed by the method of currents, which is independent 
of the parametrisation. 
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4.2 Covariant matching 

In the simplest case of field theory treated here, namely a two dimensional space- 
time, two main applications for boundary value problems (BVP) present them¬ 
selves for further discussion. Besides BVPs, initial value problems (IVP) could 
also be considered, but IVP are not of great interest for curve matching. We will 
thus forego discussing them here, although a possible application would be to 
predict the evolution of a particular model, knowing that it should “roughly” 
follow a generic model. In this case, the initial value problem must have carefully 
chosen initial values, a subject which is out of the scope of this work. The two 
applications for BVP that we will discuss are the following: 

(1) Matching between cylindrical surfaces, and 

(2) Spatio-temporal analysis. 

Case (1) In the first case, a set of slices along a cylindrical surface ( a typical 
example would be a bone) are given, where x is the parameter along the main 
axis of the surface. For the sake of simplicity, we will just consider two slices, 
but more could be added without to much trouble. The hrst step would be to 
generate the initial and final conditions, namely use un-reduction for the IVP to 
interpolate between the two curves and generate the initial and final surfaces. 
Once this is done, covariant un-reduction could be applied using a shooting 
method in time, such that the solution is a critical point of the action functional 
f L{c,Ct,Cx)dtdx. In our simple case, where the Lagrangian is purely quadratic, 
the solution would be a harmonic map, or a minimal surface, and would then 
require more advanced mathematics, beyond the present discussion. This model 
would compute the distance between two surfaces, taking into account that the 
interpolation between the slices in space should be imposed simultaneously with 
the matching in time. The resulting distance will be different than a naive model, 
which would compute the matching in time, slice by slice. For an illustration 
of matching slice by slice, we refer to the last example in [8] where a surface 
representing a nasal cavity is reconstructed out of a set of slices. The step done 
there corresponds to the generation of initial and final surfaces, whereas covariant 
un-reduction would compute the distance between these two surfaces, during the 
temporal deformation of the entire nasal cavity. 

Case (2) Spatio-temporal analysis has recently been reviewed in |1()) . from 
yet another viewpoint. Indeed, the matching in space done in [10] does not de¬ 
pend on a space parameter, but is instantaneous, namely given by a single map 
between the two curves. They also included a “time warp” which account for 
the change of pace of the evolution of the two models to be compared. In our 
case, the spatial variable comes into play on the same footing as time, and, thus, 
brings more flexibility into the comparison. Again, the theory of harmonic maps 
could help in understanding the properties of the solutions, and it is possible 
that the time warp reparameterisation could be recovered as well. 

We finally want to mention the freedom to choose the vertical force in the 
un-reduction equations. This force could be used to control the parametrisation 
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during the matching. Different types of forces could be considered, such as a 
force which would optimally redistribute the parametrisation in different regions 
of the curve, such that the number of points for discretising the curve could be 
reduced. Another force could be used to match the paramerisation of the target 
curve, such that the computationally more expensive method of currents could 
be avoided. 

Acknowledgements We are grateful to M. Bauer, M. Bruveris, L. Younes, 
S. Durrleman, R. Montgomery and T. Ratiu for valuable discussions during 
the course of this work. AA acknowledges partial support from an Imperial 
College London Roth Award, AA and DH from the European Research Council 
Advanced Grant 267382 FCCA. MCL has been partially funded by MINECO 
(Spain) under projects MTM2011-22528 and MTM2010-19111. MCL wants to 
thank Imperial College for its hospitality during the visit in which the main ideas 
of this work were developed. 

References 

1. A. Arnaudon, M. Castrillon Lopez, D.D. Holm, Covariant un-reduction and appli¬ 
cations, in preparation. 

2. M. Bauer, M. Bruveris, S. Marsland, P.W. Michor, Constructing reparameterization 
invariant metrics on spaces of plane curves, Differential Geom. Appl. 34 (2014), 
139-165. 

3. M. Bauer, M. Bruveris, P.W. Michor, Overview of the Ceometries of Shape Spaces 
and Diffeomorphism Groups, J. Math. Imaging and Vision 50 (2014), 67-90. 

4. M. Bruveris, D.C.P. Ellis, F. Gay-Balmaz, D.D. Holm, Un-reduction, Journal of 
Geometric Mechanics 3 (2011), 363-387. 

5. M. Gastrillon Lopez, T.S. Ratiu, Reduction in Principal Bundles: Covariant 
Lagrange-Poincare Equations, Comm. Math. Phys. 236 (2003), 223-250. 

6. M. Castrillon Lopez, P.L. Garcia, T.S. Ratiu, Euler-Poincare reduction on prin¬ 
cipal bundles, Lett. Math. Phys. 58 (2001), 167-180. 

7. Cendra, Hernan and Marsden, Jerrold E and Ratiu, Tudor S Lagrangian reduction 
by stages, American Mathematical Soc. 722, (2001), x-l-108 

8. C.J. Cotter, A. Clark, J. Peiro, A reparameterisation based approach to geodesic 
constrained solvers for curve matching, Int. J. Comput. Vis. 99, (2012), 103121 

9. C.J. Cotter and D.D. Holm, Geodesic boundary value problems with symmetry, J. 
Geom. Mech. 2, no. 1 (2010), 417-444. 

10. S. Durrleman, X. Pennec, A. Trouve, J. Braga, G. Gerig, N. Ayache, Toward a 
comprehensive framework for the spatiotemporal statistical analysis of longitudinal 
shape data, Int J Comput Vis. 103, no. 1 (2013), 22-59. 

11. D.C. Ellis, F. Gay-Balmaz, D.D. Holm, T.S. Ratiu, Lagrange-Poincare field equa¬ 
tions, Journal of Geometry and Physics 61, no. 11, (2011), 2120-2146. 

12. J. Glaunes, A. Trouve, L. Younes, Modeling planar shape variation via hamiltonian 
flows of curves, in H. Krim Jr. and Y.A (Eds.), Statistics and analysis of shapes 
(2006), Basel Birkhauser. 

13. D.D. Holm, J.E. Marsden and T.S. Ratiu, The Euler-Poincare equations and 
semidirect products with applications to continuum theories, Adv. in Math., 137, 
(1998), 1-81. 



Covariant un-reduction for curve matching 


13 


14. P.W. Michor, D. Mumford, An overview of the Riemannian metrics on spaces 
of curves using the Hamiltonian approach, Appl. Comput. Harmon. Anal. 23 (1) 
(2007), 74-113. 



